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' The real-time dynamics of finite-temperature gauge theories can be approximated, 

^ to leading-order accuracy in the coupling constants, by a classical field theory 

with the hard thermal loop Lagrangian. I show how this approach can be used in 
numerical lattice simulations to study dynamics of the Abelian Higgs model in or 
slightly out of equilibrium. 
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^ ■ 1 Hard thermal loops 

0^ ' Depending on its details, the electroweak phase transition may explain the 

(N ■ observed baryon asymmetry in the Universe.EI Although static equilibrium 

properties such as the phase structure of the theory have been determined to 
high accuracy with numerical lattice simulations,a much less is known about 
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' non-equilibrium dynamics, or even real-time correlators in equilibrium. The 



reason is that to calculate any real-time correlator, one needs to evaluate a 
^ • Minkowskian path integral, which cannot be done using Monte Carlo methods, 

I I and perturbation theory breaks down because of infrared problems. 

Assuming that the fields are initially in thermal equilibrium, the occupa- 
(-H I tion number of the soft modes (with momentum k<gT) is large, and they 

can be approximated by classical fieldsB The hard modes {k > T) have a small 
. ^ , occupation number and must be treated as quantum fields, but perturbation 

■ theory works well for them, and they can be integrated out perturbatively. 

H \ This is done by choosing a lattice cutoff 5x — 1/ h. with gT <ti A -^T, calcu- 

lating one-loop correlators both in the full theory and in the lattice theory and 
determining the effective Lagrangian by matching the results. Because A acts 
as an infrared cutoff, this construction is free from infrared problems, and the 
expansion parameter is g^, not g^T/ra as usually is, the case in perturbative 
calculations. In practice, since the full lattice resultu is very cumbersome and 
therefore unsuitable for numerical simulations, it is easier to match only static 
correlators, which gives a small error to the results. 

In this way, an effective theory can be constructed for the soft modes, and 
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the dynamics of this effective theory is well described by classical equations 
of motion. This approach has prjeidously been used to measure the sphaleron 
rate of hot SU(2) gauge theoryuQ I will argue that it can also be used for 
simulating non-equilibrium dynamics of phase transitions. _ 
In the case of the Abelian Higgs model, the effective Lagrangian isQ 

+ p^0)*i^^0-m|0*0-A(0»', (1) 

where mj^ = ^e^T^ - (5m|, and to| = _|_ (£2^4 ^ X/3)T'^ - Sm^. The 
mass counterterms ardj (5m|, = e^Sr/47r(5x and 5m3p = (3e^ + 4A)ST/47rfo. 
The integration is taken over the unit sphere of velocities v = (1, w), = 1. 
The extra terms in the Lagrangian depend only on the couplings and the 
temperature, as long as the high-temperature approximation T ^ m is valid. 
The equations of motion can be derived from Eq. ([l]), and are 

/fin njV.A 
-E^ - 2elmrD-cj,, 
An V ■ o 

Df,D>'(t)=-m^(l}-2X{<j)*<j))(t). (2) 

Because of the derivative in the denominator, the gauge field equation of 
motion is non-local. 



2 Local formulation 

To make numerical simulations feasible, one needs a local formulation for 
the theory. The most straightforward approaches involve describing the hard 
modes by a large number of cltarged point particles,El or by the phase-space 
distribution of these particles.BH In practice, both formulations lead to a 5-1-1- 
dimensional field theory. However, in the Abelian case, one can integrate out 
one of the dimensions, thus obtaiuing a 4-|-l-dimensional theory that is com- 
pletely equivalent with the others.0 It consists of two extra fields f(t, x, z) and 
9{t,x,z), where z e [0, 1]. In the temporal gauge, they satisfy the equations 
of motion 

d^oRz) = z^vH+ mnz^J^^V x A, 

dle{z) = z^v ■ {ye - m^i*) , 
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+mD I dzz'\Ve-mDA+^l^-^Vxfj. (3) 

With these equations, it is possible to calculate any real-time correlator at 
finite temperature. One simply takes a large number of initial configurations 
from the thermaLensemble with the probability distribution exp{—(3H), using 
the HamiltonianQ corresponding to the equations of motion (^), and evolves 
each configuration in time, measuring the correlator of interest. The average 
over the initial configurations gives the ensemble average of the correlator. 

3 Simulations 

In addition to the standard lattice discretization, the dependence on the new 
coordinate z needs to be discretized as well.Q A convenient way is to define 
the canonical momenta F = dgf and 11 = do9 and to expand both the fields 
and the momenta in terms of Legendre polynomials: 

/^"^ - dzZ^^-^P2n{z)f{z), = 1^ dzP2n{z)e{z), 

J-^P2n(z)n^), n(") = /'dzP2„Wn(z). (4) 

Jq z V t — z Jq 

The Hamiltonian can be written in terms of these Legendre modes, and 
the ensemble of initial configurations can be generated using standard Monte 
Carlo techniques. This can be done in two steps: 

(i) The Hamiltonian is Gaussian in the hard fields /, 9 and H, and they 
can therefore be integrated out analytically. This leads to the Hamilto- 
nian of the ordinary classical Abelian Higgs model, with an extra Debye 
screening term for the electric field. This classical Hamiltonian can be 
used to generate the initial configuration for the soft modes. 

(ii) Given the soft configuration from step (i) , the hard field configuration can 
be generated very effectively, since the hard Hamiltonian is Gaussian. 

After the initial configurations have been generated, any real-time correlator 
can be measured as was explained in the end of Sec. |[ 

4 Non-equilibrium dynamics 

In addition to equilibrium real-time correlators, the formulation presented 
here can also be used to study non-equilibrium dynamics, provided that the 
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hard modes remain close enough to the equiHbrium. This is plausible in 
a phase transition, since the phase of a system is a property of the long- 
wavelength modes only. Near the phase transition, in both phases, all the 
masses are suppressed by powers of the coupling constant g relative to the 
temperature, and the high-temperature approximation can be used. In this 
approximation, the phase of the system does not enter the results of the one- 
loop diagrams, i.e. the distribution of the hard modes is indeed the same in 
both phases, within the accuracy of our approach. Thus, there is no reason 
for the hard modes to fall out of equilibrium during the transition. 

A natural way of studying a phase transition would be to start from ther- 
mal equilibrium in the Coulomb phase, and decrease the temperature so that 
the system undergoes a transition to the Higgs phase. This requires a mech- 
anism for changing the temperature, and in practice, it is easier to keep the 
temperature constant and change the parameters, such as the mass of the 
Higgs field, instead. In fact, when A <^ e^, even that is not necessary. The 
transition is of first order, and if one thermalizes the system initially to the 
metastable Coulomb phase below Tc, bubbles of the Higgs phase nucleate dur- 
ing the time evolution, and the phase transition takes place. Assuming that 
the latent heat is small enough, the temperature does not change significantly. 
In this way, many non-equilibrium properties of the phase transition can be 
studied non-perturbatively. 
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